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ABSTRACT 

Theoretical models of very metal-poor intermediate-mass Asymptotic Giant Branch (AGB) stars predict a 
large overabundance of primary nitrogen. The very metal-poor, carbon-enhanced, s-process-rich stars, which 
are thought to be the polluted companions of now-extinct AGB stars, provide direct tests of the predictions 
of these models. Recent studies of the carbon and nitrogen abundances in metal-poor stars have focused on 
the most carbon-rich stars, leading to a potential selection bias against stars that have been polluted by AGB 
stars that produced large amounts of nitrogen, and hence have small [C/N] ratios. We call these stars Nitrogen- 
Enhanced Metal-Poor (NEMP) stai's, and define them as having [N/Fe] > +0.5 and [C/N] < -0.5. In this 
paper, we report on the [C/N] abundances of a sample of 21 carbon-enhanced stars, all but three of which 
have [C/Fe] < +2.0. If NEMP stars were made as easily as Carbon-Enhanced Metal-Poor (CEMP) stars, then 
we expected to find between two and seven NEMP stars. Instead, we found no NEMP stars in our sample. 
Therefore, this observational bias is not an important contributor to the apparent dearth of N-rich stars. Our 
[C/N] values are in the same range as values reported previously in the literature (-0.5 to +2.0), and all stars 
are in disagreement with the predicted [C/N] ratios for both low-mass and high-mass AGB stars. We suggest 
that the decrease in [C/N] from the low-mass AGB models is due to enhanced extra-mixing, while the lack of 
NEMP stars may be caused by unfavorable mass ratios in binaries or the difficulty of mass transfer in binary 
systems with large mass ratios. 

Subject headings: nuclear reactions, nucleosynthesis, abundances-stars: abundances-stars: AGB and post- 
AGB-stars:carbon-stars: Population II 



1. INTRODUCTION 

Very Metal-Poor (VMP) stars ([Fe/H] < -2.0)^ provide es- 
sential tools for the study of element production in stars and 
galactic chemical processing during the early stages of the 
evolution of our Galaxy. The number of VMP stars with mea- 
sured elemental abundance ratios has been increasing rapidly 
in the past decade. The large modern surveys for metal-poor 
stars, most importantly the HK survey of Beers and colleagues 
(Beers et al. 1992; Beers 1999) and the Hamburg/ESO Sur- 
vey (HES) of Christlieb and collaborators (Christlieb 2003), 
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have produced medium-resolution confirmation spectra of 
many thousands of metal-poor candidates; the most interest- 
ing of these have been (and are being) followed up with high- 
resolution spectroscopic studies on large-aperture telescopes 
(e.g. Cayrel et al. 2004; Honda et al. 2004a; Johnson & Bolte 
2004; Barklem et al. 2005). 

One of the most important discoveries of these new sur- 
veys is that at least 20 % of all VMP stars exhibit conspicu- 
ous enrichments of the CNO elements, most notably C (Beers 
& ChristHeb 2005; Lucatello et al. 2006). These Carbon- 
Enhanced Metal-Poor (CEMP) stars are defined by Beers & 
Christlieb as metal-poor stars with [C/Fe] > +1.0. Many of 
the CEMP stars exhibit overabundances of the elements as- 
sociated with s-process nucleosynthesis. Aoki et al. (2003) 
estimate this fraction to be between 70 % and 80 %; Beers & 
Christlieb refer to these as the CEMP-s stars. The observed 
abundance patterns for CEMP-s stars suggest nucleosynthetic 
origin in low- or intermediate-mass stars that have evolved 
through the thermally pulsing asymptotic giant branch (AGB) 
phase, and later transferred this processed material to a surviv- 
ing low(er)-mass companion. In this sense the CEMP-s stars 
are the metal-poor analogs of the classical CH stars (Keenan 
1942). The AGB star that originally enriched the presently 
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Fig. 1 . — Surface abundances of AGB stars of different masses for metal- 
licity [Fe/H] = —2.3 over several dredge-up episodes. As more dredge-up 
occurs, the surface abundance of C increases. For the higher-mass cases, 
dredged-up C is processed into N via HBB, leading to lower C/N ratios. The 
final total yields of the elements produced by the AGB stars and ejected into 
the ISM by mass loss are given in Herwig (2004b). 

observed companion VMP star is now a white dwarf; its pres- 
ence is often revealed by tell-tale radial-velocity variations 
observed for the companion star. Based on the still-limited 
numbers of CEMP-s stars for which sufficient multi-epoch 
spectroscopic data has been obtained, a 100% binary frac- 
tion is possible (Lucatello et al. 2005). It is thus expected that 
the large overabundances of a number of elements observed in 
CEMP-s stars should reflect the AGB nucleosynthetic yields, 
perhaps with some modification due to giant-branch evolution 
of the companion. 

In all AGB stars C is produced by the triple-a reaction in 
He-shell flashes and convectively dredged up into the stellar 
envelope. For stars of low initial mass (2-3 M0), this even- 
tually leads to a C-rich composition, with C/O > 1. In these 
models N is not enhanced because it is burned during the He- 
shell flash. During the interpulse phase the convective enve- 
lope has no mixing connection with the H-shell, according 
to standard models, and no alteration of N is expected. For 
larger masses (> 3.5 M0, depending on metallicity) the effi- 
ciently dredged up carbon is transformed into N by the hot- 
bottom burning (HBB) process (Lattanzio 1992; Boothroyd 
et al. 1993; Forestini & Charbonnel 1997; Herwig 2004a; 
Ventura & D'Antona 2005). As a result, the models predict 
small C/N ratios in the stellar envelope ([C/N] « -1). 

Simulations of low- and intermediate-mass stellar evolution 
confirm that the well-established trends of C and N produc- 
tion in solar-metallicity AGB stars extend to very low metal- 
licity (Herwig 2004b). The amount of N produced is indepen- 
dent of the initial metallicity of the star, because it is based on 
the primary production of C in the He-burning shell. There- 
fore, HBB in intermediate-mass VMP AGB stars provides a 
primary source of nitrogen in the early Universe. Figure 1 
shows the evolution of [C/N] ratios on the surface of metal- 
poor AGB stars from the models of Herwig (2004b). Two ef- 
fects are immediately clear: (1) a lower predicted [C/N] ratio 
in more massive AGB stars, and (2) a lower predicted [C/Fe] 
on the surfaces of these same stars. 

The predicted time-averaged yields for AGB stars are 
shown in Figure 2. In this Figure we also plot [C/N] ratios 
for CEMP stars that were available in the literature when we 
began this project. We would expect that at least one-third 
of the stars should have been polluted by intermediate-mass 



AGB stars, assuming (1) stars between 3.5 Mq and 7.5 Mq 
undergo HBB, (2) stars between 1 .0 M0 and 7.5 Mq produce 
enough C to enhance their companions during mass transfer, 
and (3) the binary mass ratios from Pinsonneault & Stanek 
(2006) are correct for low-metallicity stars. It is evident that 
the observed [C/N] ratios in CEMP stars do not fit the pre- 
dictions of either the low-mass or the intermediate-mass AGB 
models, but instead fall in an intermediate regime that is not 
covered by the models. 

We are left with two questions. First, why do the observa- 
tions of stars with large [C/Fe] show larger N (smaller [C/N]) 
than predicted by models for the evolution of 2-3 Mq stars? 
This problem is reminiscent of the high N abundances ob- 
served in the moderately metal-poor CH stars, and may be 
related to mixing processes that are not accounted for in mod- 
els of AGB evolution (Vanture 1992). This issue will be 
discussed later in this paper. Second, where are the VMP 
stars that were, in the context of the mass-transfer scenario, 
polluted by the N-rich 3.5-7.5 Mq stars? The solution for 
this problem is less obvious. The HBB responsible for the 
very efficient production of N in intermediate-mass stars is 
a robust prediction of stellar-evolution models. If the binary 
mass-transfer scenario is invoked for the CEMP-s stars with 
[C/N] sa 1, then one may ask why there are no NEMP stars 
with [C/N] « -1, as might be expected to arise from systems 
in which the donor is a more massive HBB AGB star. Before 
considering possible solutions to this second problem we need 
to test one obvious possibility, that the absence of NEMP stars 
among previously analysed VMP stars may be simply the re- 
sult of a selection bias against their detection. 

Because of the requirement for large-aperture telescopes 
in order to obtain high-resolution, high-S/N spectra of the 
CEMP stars in the HK and HES, only a limited number of 
these objects have been observed to date. Special attention 
has been focused on the subsample of CEMP stars that are 
likely to be the most iron-poor and/or carbon-rich (e.g. Nor- 
ris et al. 1997; Aoki et al. 2002d). This may have led to an 
observational bias against the discovery of NEMP stars. The 
CH features at 4305 A are routinely covered in the medium- 
resolution spectroscopic follow-up of metal-poor candidates 
from the HK survey and the HES; thus, CEMP stars are read- 
ily identified and placed on target lists for examination at high 
resolution. However, N-rich stars are not so easily recog- 
nized. The medium-resolution confirmation spectra generally 
extend no bluer than about 3600 A; the only N-sensitive fea- 
ture included in this range is the CN band at 3850 A, which 
requires large enhancements of both C and N to be strong. As 
a result, the majority of stars in the current literature with de- 
tailed studies of their elemental abundance patterns are those 
with [C/N] > 0, i.e., they are more C-rich than they are N- 
rich. Thus, the lack of stars with [C/N] < could well be 
an observational selection effect, especially since the lower 
C-enhancements predicted in intermediate-mass AGB stars 
would make their companions potentially less C-rich, and 
hence less likely to be followed up at high spectral-resolution. 

Previous low -resolution surveys of metal -poor ([Fe/H] < 
-1.0) dwarf stars that included the blue NH band have dis- 
covered some stars with stronger than average nitrogen abun- 
dances, including HD 74000 and HD 160617 (Bessell & 
Norris 1982), HD 97916 and HD 166913 (Laird 1985) and 
HD 25329 (Carbon et al. 1987). Higher-resolution foflow- 
up studies (e.g., Beveridge & Sneden 1994; Mashonkina et al. 
2003) have shown that these stars are somewhat s-process rich 
([Ba/Fe] ~+0.5). They do not appear to be C-enhanced, how- 
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Fig. 2.— Literature values for [C/N] ratios in VMP as of 2002. The lines 
are labeled with the masses of the AGB stars from the models of Herwig 
(2004b). We assume that the initial C and N on the observed star's suiface 
are overwhelmed by the contribution from the AGB star However, because 
the [C/Fe] ratio depends on both the amount of carbon transfeiTed from the 
AGB star and the amount of Fe on the companion star surface, there is no a 
priori dependence of the [C/Fe] ratio on the mass of the AGB star. We do 
not include in this figure CS 22892-052, which is r-process rich, CS 22949- 
037, which is extremely O-rich, orCS 30314-067, which has new abundances 
reported in this paper. 

ever, and their binary status is not confirmed. These stars are 
probably related to the phenomenon discussed here, but we 
confine ourselves in the rest of the paper to the discussion of 
CEMP stars found in recent surveys. 

We have undertaken a medium-resolution observing cam- 
paign to address the question of whether observational selec- 
tion biases might be responsible for the lack of known NEMP 
stars. We examine a sample of VMP stars with moderate 
carbon enhancements, -1-0.5 < [C/Fe] < -l-l.O (based on their 
medium-resolution confirmation spectra), in order to better 
constrain the range of the N enhancements in these stars, and 
to see if the N abundances obtained agree better with expec- 
tations based on AGB models. This is accomplished using 
near-UV medium-resolution spectroscopy that covers the re- 
gion of the NH band at 3360-3370 A for a sample of 21 mod- 
erately carbon-enhanced metal-poor stars, along with a num- 
ber of similar stars from the literature with available high- 
resolution results. In § 2 we describe our sample selection 
criteria, observations, and data reduction procedures. Details 
of our abundance analysis for this sample are provided in § 3. 
In § 4 we summarize our results. A discussion of the theoret- 
ical expectations and a comparison with our present results is 
provided in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection 

We wished to observe a sample of stars with [Fe/H] < -2.5 
and [C/Fe] between -1-0.5 and -l-l.O. Our primary source of 
targets was the list of [Fe/H] and [C/Fe] values for the HK 
survey stars calibrated by Rossi et al. (2005). We also added 
stars from the Beers, Preston, & Shectman (1992) list of stars 
with abnormally strong G-bands (their Table 8). To com- 
pare our derived abundances with those obtained from higher- 



resolution spectra, we also took spectra of some well-known 
bright metal-poor HD and BD stars, as well as three stars from 
the HK survey, namely CS 22892-052, CS 22968-014, and 
CS 22950-46. These three are not part of our C-enhanced 
sample either, because they are known to be r-process-rich 
and probably polluted by a different mechanism altogether 
(CS 22892-052) or are not C-enhanced (CS 22968-014 and 
CS 22950-048). Ideally, all of our stars would be subgiant or 
main-sequence stars, to avoid possible alterations in the sur- 
face C and N abundances due to CN processing and mixing on 
the red giant branch, which can result in decreasing C and in- 
creasing N abundances (Gratton et al. 2000; Spite et al. 2005). 
However, molecular features such as CH and NH are much 
weaker in hotter stars, so we decided to include giants in our 
list. Any CN cycling that might have operated would serve 
to increase the N abundances, and therefore bias the study to- 
ward finding more NEMP stars than would otherwise be the 
case. This turned out not to be a concern in the final analysis. 

2.2. Observations 

To observe the NH band we require near-UV sensitivity, 
however, a resolution of 2-2.5 A is more than adequate. The 
RC spectrographs at Kitt Peak National Observatory and at 
Cerro Tololo Interamerican Observatory have these capabili- 
ties, allowing us to observe both northern and southern hemi- 
sphere targets. We observed 8 stars at KPNO over three 
nights from 30 Aug to 1 Sep 2003. We employed the BL420 
grating and the CUSO4 order blocking filter. This produced 
a wavelength coverage of 3250 A to 4700 A with a resolu- 
tion of 2. OA. The F3B chip was used, with a gain setting 
of 2.3 e"/ADU, and a read noise of 7.5 e". The dispersion 
was 0.76A/pix. A quartz lamp with a BG3 blocking filter 
was used to obtain flatfields, and the FeAr lamp was used 
for wavelength calibration. We also observed 20 stars over 
6 nights at CTIO in two observing sessions: 6-8 Sep 2003 
and 26-28 Nov 2004. Some part of each run was lost to 
weather We used the RC Spectrograph with the blue collima- 
tor, the KPGLl grating, and the CUSO4 order-blocking filter 
We adopted Decker 2, which provided a 1 arcsec slit. This 
produced a wavelength coverage of 3250 A to 4700 A with a 
resolution of 2.5 A. The gain setting was 1.94e~/ADU, and 
the read noise was 7.5 e". The dispersion was 0.95A/pix. A 
helium-neon-argon lamp was used for wavelength calibration 
for these data. The second run at CTIO (26-28 Nov 2004) 
was plagued by large amounts of scattered light in the spec- 
tra. Although the scattered light was subtracted, the Poisson 
noise from its presence reduced the S/N around the NH band 
to ^ 10-20. We were still able to obtain useful upper limits on 
N, and therefore lower limits on [C/N], with these data. Ta- 
ble 1 present a summary of the observations. Figure 3 shows 
a typical spectrum obtained from these observations. 

2.3. Data Reduction 

The data were reduced using IRAF^. The usual reduction 
procedures were applied: bias-subtraction, flat-fielding, spec- 
tral extraction, sky subtraction, and wavelength calibration. 
For most stars, with the exception of the observations in 2004, 
we achieved our desired S/N of 50/1 per resolution element. 
The S/N at the NH bandhead is listed in Table 1. 

IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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TABLE 1 
Log of Observations 



Star 


V 


Telescope 


Date of Obs. 


Exposures 


S//Vat3400A 


BD-18°5550 


9.270 


CTIO 


2003 Sep 05 


2x600 


192 


HD 122563 


6.196 


KPNO 


2003 Aug 31 


3 X 200s 


267 


HD 160617 


8.740 


CTIO 


2003 Sep 04 


3 X 600s 


207 


HD 186478 


8.920 


CTIO 


2003 Sep 04 


3 X 600s 


425 


CS 22174-007 


12.409 


CTIO 


2003 Aug 30 


3 X 1800s 


142 


CS 22183-031 


13.622 


CTIO 


2004 Nov 28 


3 X 1200s 


17 


CS 22879-029 


14.425 


CTIO 


2003 Sep 06 


Ixl200s,2x900s 


117 


CS 22884-097 


14.868 


CTIO 


2003 Sep 05 


3 X 1200s 


33 


CS 22887-048 


12.866 


KPNO 


2003 Aug 30 


Ixl800s,2x900s 


123 


CS 22891-171 


14.293 


CTIO 


2003 Sep 06 


3 X 1200s 


58 


CS 22892-052 


13.213 


CTIO 


2003 Sep 04 


4 X 600s 


115 


CS 22898-062 


13.788 


KPNO 


2003 Aug 31 


3 X 1800s 


71 


CS 22945-024 


14.360 


CTIO 


2003 Sep 04 


3 X 1200s 


55 


CS 22947-187 


12.962 


CTIO 


2003 Sep 04 


3 X 1200s 


117 


CS 22948-104 


13.929 


CTIO 


2004 Nov 28 


3 X 1200s 


17 


CS 22949-008 


14.168 


KPNO 


2003 Aug 31 


1 X 600s 3 X 1800s 


95 


CS 22950-046 


14.224 


KPNO 


2003 Aug 30 


4x 1800s 


45 


CS 22958-042 


14.516 


CTIO 


2003 Sep 06 


4 X 1200s 


42 


CS 22958-083 


14.423 


CTIO 


2044 Nov 27 


3x1200 


17 


CS 22960-053 


14.830 


CTIO 


2003 Sep 06 


3 X 1200s 


40 


CS 22968-014 


13.684 


CTIO 


2003 Sep 04 


3 X 1200s 


80 


CS 29493-090 


14.039 


CTIO 


2003 Sep 04 


3 X 1200s 


57 


CS 29495-042 


14.516 


CTIO 


2003 Sep 06 


1 X 1200s 2 X 900s 


67 


CS 29497-030 


12.656 


CTIO 


2003 Sep 06 


Ixl200s,3x600s 


134 


CS 29512-073 


14.137 


KPNO 


2003 Aug 30 


2xl800s,lx900 


66 


CS 30314-067 


11.817 


CTIO 


2003 Sep 04 


3 X 1200s 


122 


CS 31062-041 


13.934 


KPNO 


2003 Aug 31 


3 X 1800s 


53 


CS 31080-095 


12.989 


CTIO 


2004 Nov 27 


4 X 1200s 


17 



3.1. Molecular and Atomic Data 

We assembled lineUsts covering 3 340-3400 A for NH and 
4270-4330 A for CH. The atomic parameters are from the Vi- 
enna Line Database (Piskunov et al. 2000). The molecular 
line lists were in general adopted from Kurucz (2006b) and 
are discussed further below. We adopted the solar values of 
Anders & Grevesse (1989), with the exception of Fe, where 
we assume a solar value of log e = 7.52. While the solar values 
of C and N have been the subject of much dispute in recent 
years (Asplund et al. 2005), in our case they merely represent 
scahng factors apphed to both the yields from AGB stars of 
Herwig (2004b), and to our derived loge for C and N from 
these stars. 

3.1.1. NH 

We used the bandhead for the A-X (0-0) and (1-1) transi- 
tion at 3360 A and 3370 A, respectively. The only important 
isotope is '^N, because the '^N/'^N ratio is > 100 wherever 
it has been measured (Chin et al. 1999). The wavelengths 
of NH from Kurucz (2006b) are in excellent agreement with 
the laboratory wavelengths measured by Brazier et al. (1986) 
for the regions near the NH bandhead. However, Shavrina 
et al. (1996) pointed out that the Kurucz gf- values were too 
high by a factor of two. Therefore, we have adopted Kurucz 
wavelengths, but divide his oscillator strengths by this factor. 
These gf-values are then on the same scale as those of Sneden 
(1973). 

Conflicting values for the dissociation potential of NH ex- 
ist in the hterature. Seal & Gaydon (1966) measured 3.21 ± 
0.16eV. More recent experiments have found higher values. 
Graham & Lew (1978) measured Do < 3.47 ±0.05 eV, while 
Tarroni et al. (1997) measured Do > 3.419 iO.OlOeV. With 
these results in mind, we have adopted 3.45 eV, close to the 
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Fig. 3.— A normalized spectrum of CS 22947-187 from the RCSpec at 
CTIO. 



3. ABUNDANCE ANALYSIS 

We used the 2002 version of MOOG (Sneden 1973) for our 
synthesis of the NH, CH and Ca II K regions. We interpolated 
the grids of Kurucz (2006a), using the models with overshoot- 
ing, but without a-enhancements or the new opacity distribu- 
tion functions, because that grid was the most extensive at the 
lowest metallicities. For the accuracy we require in [C/N] at 
these low metalUcities, our choice of model atmosphere grid 
is not a significant source of error. 
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Fig. 4. — Synthesis of the NH bandhead region in HD 186478 for a high- 
resolution spectrum (top) and a low-resolution spectrum (bottom). The model 
atmosphere used in this paper was adopted for both syntheses. One line 
shows the synthesis without any N, while the second shows the synthesis 
with loge(N)=6. 12. Filled squares show the data. 

value used by Spite et al. (2005) in their study of C, N and O 
abundances in VMP stars. 

Our linelist in the NH region is illustrated in Figure 4. We 
first show the fit to a high-resolution spectrum of the metal- 
poor star HD 186478 from Johnson (2002) and then to the 
lower-resolution spectrum taken for this study. 

If we used a dissociation potential of 3.21 eV instead of 
3.45 eV, our derived N abundances would increase by ^ 
0.3 dex. If we used the uncorrected Kurucz gf-values, our 
abundances would decrease by a similar amount. The lack of 
agreement among different studies on the dissociation poten- 
tial and oscillator strengths is a major source of the differences 
in the derived C and N abundances in the present literature. 
Fortunately, these issues are dwarfed by the > 1 dex differ- 
ence we hope to observe. 

3.1.2. CH 

Our analysis of the carbon abundances for program stars 
used the CH G band feature. Again, we used the Kurucz line 
list for this region. The wavelengths and oscillator strengths 
agreed well with the more limited lists of wavelengths from 
Zachwieja (1995, 1997) and oscillator strengths from LIF- 
BASE (Luque & Crosley 1996). Theoretical models predict 
different '^C/^^ ratios for AGB stars of different masses, rang- 
ing from > 1000 for 2Mp, to - 5 for 6M0. The G band 
contains both '^CH and CH lines. However in our low- 
resolution spectra, these lines are intermingled to such an ex- 
tent that if the total carbon abundance does not change, our 
synthesis of the G-band remains essentially independent of 
the carbon isotope ratios. We adopted a value of 80 for our 
synthesis. Figure 5 shows the fit in the G-band region for 
HD 186478. 

3.1.3. C and N Abundances from Different Molecular 
Species 

Most studies of metal-poor stars measure the N abun- 
dance from the blue CN system near 3870-3880 A, as this 
wavelength region is usually covered in the setups used for 
medium-resolution spectroscopy. Spite et al. (2005) were able 
to work with high-resolution spectra of metal-poor stars that 
covered both the NH and CN features. They found that the 
NH lines consistently yielded 0.3 dex higher N abundances 
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Fig. 5. — Synthesis of the CH bandhead region in HD 186478 for a high- 
resolution spectrum (top) and a low-resolution spectrum (bottom). The model 
atmosphere used in this paper was adopted for both syntheses. One line 
shows the synthesis without any C, while the second shows the synthesis 
with loge(C) = 5.73. Filled squares show the data. 

than those obtained from the CN lines. The source of the dis- 
agreement is unclear (and we do not resolve it here), but note 
that our N abundances (derived from NH) would be lower by 
0.3 dex if the CN abundance scale were the correct one. Lit- 
erature sources also derive C abundances from both the C2 
and the CH features. Aoki et al. (2002a) measured a higher C 
abundance, by 0.2 dex, in the subgiant LP 625-44 from the C2 
lines as compared to the CH feature. These differences could 
result in systematic offsets in C or N abundances; however, 
their magnitude is smaller than the size of the expected [C/N] 
range. Where necessary for our discussion, we thus adopt lit- 
erature C and N values as stated, selecting C from CH and N 
from NH whenever possible. 

3.2. Model Atmosphere Parameters 

We used broadband photometry, provided by either Beers 
et al. (2006) or taken from the SIMBAD listing, to derive Teff 
and isochrones to derive gravities. We also took advantage of 
the well-known correlation between log g and ^ for estima- 
tion of microturbulent velocities, as described below. Finally, 
we synthesized the Ca II K line to determine the appropriate 
metallicity. We discuss additional details and comparisons 
with other methods of deriving model atmosphere parameters 
below. 

3.2.1. Teff 

We used the color-Teff relations from Alonso et al. (1996, 
1999). Accurate photometry for our targets in Johnson UBV, 
and Kron-Cousins RI was obtained as part of an ongoing pro- 
gram of observing metal-poor stars (Beers et al. 2006). Near- 
IR JHK magnitudes were retrieved from the 2MASS catalog 
(Skrutskie et al. 2006). Table 2 lists the available photome- 
try. We have adopted reddening values from Schlegel et al. 
(1998). Almost all our targets have distances > 1 kpc, much 
larger than the llOpc scale height of the dust (Mendez & 
van Altena 1998), hence the assumption that all the redden- 
ing is between us and the star is appropriate. For the more 
nearby stars, we calculated the correct amount of reddening 
to use based on the model of Mendez & van Altena (1998) 
(their Equation 4), iterating until the following quantities con- 
verged: assumed reddening based on distance from Sun, Teff 
from V -K, log g and My from isochrones, and finally dis- 
tance from My. This required between one and three cycles. 
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TABLE 2 
Photometry Data 



Star 


V 


B-V 


U-B 


V-R 


V-I 


V-J 


V-K 


E(B-V) 




9 270 


9 270 


910 


{) 2-^{) 


530 


2 068 


2 715 


17 




6 196 


912 


0.340 


0.580 




1.815 


2 473 


02 




O. I^\J 










1112 


1 429 


01 


n.i-' i OUT- / o 


8 920 


930 








1.803 


2 477 


10 


79174-007 


12 409 


694 




448 


900 


1 524 


2 091 


03 


V.I Z.Z.iOJ WJl 


13 622 


673 




444 


898 


1511 


2 040 


04 


79SS4 nQ7 


14 868 


579 




0.407 


849 


1 401 


1 742 


23 




1 7 Sft6 


u. joy 




97S 


556 


n QOi 


1 192 


OS 


PS 79S7Q 070 


14 425 






980 


604 


1 000 


1 9QQ 


04 


PS 79S^Q1-171 


14 293 


0.872 




485 


918 


1.522 


2.030 


0.07 


PS 79SQ7-0'^9 

V.I Z.Z.07i. WJi. 


13 213 


800 




491 


990 


1 721 


2 284 


03 


PS 79SQS-0fi9 


13 788 


627 




420 


840 


1 439 


1 937 


05 


P<; 79Q4S 094 


If . JOU 


n 799 










1 Q9^ 


0'^ 


PC 99Q47_1 R7 


12.962 


0.648 




0.454 


0.877 


1.494 


1.934 


0.09 


PC 99048-104 


13.929 


0.604 




0.410 


0.821 


1.431 


1.917 


0.02 


pc 99040-OOS 


14.168 


0.494 




0.310 


0.623 


1.045 


1.380 


0.04 


pc 990^0-046 


14.224 






0.585 


1.158 


1.974 


2.641 


0.06 


CS 22958-042 


14.516 


0.479 




0.294 


0.614 


1.001 


1.303 


0.02 


CS 22958-083 


14.423 


0.664 


-0.030 


0.442 


0.902 


1.477 


2.029 


0.04 


CS 22960-053 


14.830 


0.760 


0.100 






1.532 


2.069 


0.01 


CS 22968-014 


13.684 


0.755 




0.472 


0.969 


1.658 


2.217 


0.01 


CS 29493-090 


14.039 


0.835 




0.516 


1.028 


1.740 


2.366 


0.03 


CS 29495-042 


13.607 


0.626 




0.415 


0.825 


1.261 


1.737 


0.03 


CS 29497-030 


12.656 


0.299 




0.215 


0.440 


0.694 


0.911 


0.02 


CS 29527-048 


14.833 


0.455 




0.314 


0.632 


1.046 


1.356 


0.02 


CS 29512-073 


14.137 


0.568 




0.377 


0.753 


1.235 


1.624 


0.05 


CS 30314-067 


11.817 


1.123 




0.642 


1.234 


2.058 


2.789 


0.07 


CS 31062-041 


13.934 


0.820 




0.458 


0.950 


1.545 


2.113 


0.03 


CS 31080-095 


12.989 


0.521 


-0.291 


0.317 


0.617 


1.050 


1.394 


0.01 



We adopt the V-K temperature because it is independent 
of metalHcity and not strongly affected by carbon molecular 
features, contrary to the way that bluer colors, such as, B-V, 
have been shown to behave (e.g. Preston & Sneden 2001; 
Cohen et al. 2006). A comparison of T^aiV -K) with tem- 
peratures derived from other colors shows an average offset 
of Teff(V - K) < Teff(J - K) by 53 K, > Teff(B - V) by 60 K, 
> Teff(U- V) by 95K and > Teff(V-I) by 74 K. The higher 
temperatures for V-K compared to B-V and fZ-V are ex- 
pected if carbon absorption in the blue bands makes the stars 
appear redder than the non C-rich calibration sample. A com- 
parison with literature values for our sample (Table 4) shows 
that our Teff estimates are hotter, on average, by 60 K, with an 
rms scatter of 150 K. A reasonable random error for our Tgff 
determinations is therefore 150 K. 

3.2.2. logg 

We adopted the Z = 0.0001 ([Fe/H] = -2.31, [a/Fe] = 
-1-0.30), 14Gyr isochrone from Bergbusch & VandenBerg 
(2001). Note that the age is unimportant as long as it is greater 
than ~ 10 Gyr. In choosing a log g for a given Tgff, we assume 
that the star is a turnoff, subgiant, or giant star, as appropriate 
for its location on the isochrone. In general, this assumption is 
a good one. However, a few stars could be horizontal-branch 
stars, and, indeed, a high-resolution analysis by Mc William 
et al. (1995a) showed that CS 22947-187has the lower gravity 
of a red horizontal-branch star, logg = 1.30, compared to our 
value of logg = 3.44. We recalculated the abundances using 
the logg and ^ value ofMcWilhametal. (1995a). The derived 
[Ca/H] did not change, while the loge(N) and the loge(C) in- 
creased by 0.9 dex and 0.80 dex, respectively. The [C/N] ra- 
tio, however, changed by a mere 0.10 dex. Thus, for the small 
subset of stars which we may mistakenly classify as subgiants 



rather than horizontal-branch stars, the change in [C/N], the 
crucial aspect of this paper, is smaller than that which could 
arise from other sources of error. 

The major uncertainty in estimation of log g for almost all 
stars is the uncertainty in Teff, because a change of ±150K 
leads to a change of ±0.14 dex in log g for a turnoff star, 
0.19 dex for a subgiant and 0.35 dex for a giant. A compar- 
ison with literature values (Table 4) yields an average differ- 
ence of 0.20 dex with an rms scatter of 0.23 dex. Therefore, 
we adopt a log g error of 0.3 dex for turnoff and subgiant stars 
and 0.4 dex for giants, respectively. 

3.2.3. Microturbulent velocity 

Our abundances derived from lower dispersion data are not 
very sensitive to ^, the microturbulent velocity, but it still 
needs to be known to within 1 kms~' . Fortunately, there exists 
a well-known correlation between log g and ^. The HERBS 
collaboration (Barklem et al. 2005) measured ^ for 254 metal- 
poor stars from high-resolution spectra. We fit their log g and 
^ values with a second-order polynomial and use the relation 
for our sample: 

C = 2.822 -0.669 logg -I- 0.080 (logg)2 kms"' (1) 

This equation is valid for the range of surface gravities ex- 
hibited by stars in the HERBS sample (1.0 < logg < 4.2), 
which encompasses the expected range of surface gravity for 
our program stars. The HERBS data exhibit an rms scat- 
ter of 0.17kms"' around this relation. We tested changes of 
±0.3 km s"^ and found no change in our synthesis; therefore 
uncertainties in ^ do not contribute significantly to our final 
abundance errors. 

3.2.4. [Fe/H] 



TABLE 3 
Model Atmosphere Parameters 



Star Eg 


-V (mag) 


Teff (K) 


logg 


[m/H] 


5(km/s) 


BD-18° 5550 


0.17 


4806 


1.72 


-2.89 


1.91 


HD 122653 


0.02 


4615 


1.27 


-2.47 


2.10 


HD 160617 


0.01 


5882 


3.69 


-1.69 


1.45 


HD 186478 


0.10 


4831 


1.78 


-2.63 


1.89 


CS 22174-007 


0.03 


5059 


2.37 


-2.00 


1.69 


CS 22183-031 


0.04 


5196 


2.76 


-2.79 


1.59 


CS 22879-029 


0.04 


6300 


3.90 


-1.93 


1.43 


CS 22884-097 


0.23 


6460 


4.00 


-1.94 


1.43 


CS 22887-048 


0.05 


6455 


3.99 


-2.22 


1.43 


CS 22891-171 


0.07 


5297 


3.07 


-2.45 


1.52 


CS 22892-052 


0.03 


4861 


1.86 


-2.76 


1.86 


CS 22898-062 


0.05 


5309 


3.10 


-1.74 


1.52 


CS 22945-024 


0.03 


5289 


3.04 


-2.26 


1.53 


CS 22947-187 


0.09 


5489 


3.44 


-2.25 


1.47 


CS 22948-104 


0.02 


5270 


2.99 


-2.39 


1.54 


CS 22949-008 


0.04 


6144 


3.82 


-1.92 


1.44 


CS 22950-046 


0.06 


4604 


1.25 


-3.29 


2.11 


CS 22958-042 


0.02 


6224 


3.86 


-2.65 


1.44 


CS 22958-083 


0.04 


5189 


2.74 


-2.50 


1.59 


CS 22960-053 


0.01 


5061 


2.38 


-3.08 


1.68 


CS 22968-014 


0.01 


4892 


1.93 


-3.30 


1.83 


CS 29493-090 


0.03 


4739 


1.56 


-2.82 


1.97 


CS 29495-042 


0.03 


5400 


3.32 


-2.30 


1.49 


CS 29497-030 


0.02 


7163 


4.20 


-2.20 


1.43 


CS 29512-073 


0.05 


5751 


3.62 


-2.10 


1.45 


CS 30314-067 


0.07 


4476 


0.96 


-2.67 


2.25 


CS 31062-041 


0.03 


5042 


2.32 


-2.30 


1.70 


CS 31080-095 


0.01 


5972 


3.73 


-2.75 


1.44 



We require an estimate of [Fe/H] for two reasons, most im- 
portantly to have the correct metallicity atmospheric model, 
and secondly in order to measure the [C/Fe] ratio associated 
with the overall enhancement from AGB-star pollution. For 
these purposes, an accuracy of 0.3 dex is sufficient. Test syn- 
theses of the NH and CH features and the Ca II K lines re- 
vealed only a small (< 0.05 dex) dependence of these quanti- 
ties on the metallicity of the model atmosphere. Furthermore, 
the carbon-enhanced nature of these stars has already been 
established by larger-scale surveys with many stars to pro- 
vide a control sample. The lack of dependence on the overall 
metallicity of the model is fortunate, as we did not use carbon- 
enhanced models, in keeping with most of the previous work 
in the field. Masseron (2005, priv. comm.) has shown that 
errors on the abundances from neglecting C-enhancement is 
mitigated in the hotter, less C-rich stars that comprise our 
present sample. 

We first measured [Ca/H] using the Call K line that is 
prominent in all of our medium-resolution spectra, regardless 
of temperature and metallicity. The linelist for this exercise 
was taken from Castelli & Kurucz (2003). We then converted 
to [Fe/H] using an assumed [Ca/Fe] ratio of +0.30 dex, the av- 
erage [Ca/Fe] ratio determined by Cayrel et al. (2004) in their 
sample of non carbon-enhanced metal-poor stars. 

Several of our targets have been observed at high spectral 
resolution by previous authors. These include some of the 
most C-rich stars in our sample, as well as three bright metal- 
poor stars (HD 122563, HD 186478 and BD-18°5550) that 
we observed in order to compare with previous results. In 
Table 4, we compare our model atmosphere parameters with 
values taken from the literature, restricting ourselves to recent 
papers in the case of the HD and BD stars. In general, our 
model atmospheres have somewhat higher Teff and higher log 
g values. Our [Fe/H] values are also higher, which is expected 



I 0-2 ^ ^ 
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lambda (A) 
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^ ' ' ' 

4300 4290 4300 4310 

lambda (A) 

Fig. 6.— Example of syntheses for CS 22891-171 for (top) NH and (bot- 
tom) CH. In each case, the solid line represents the adopted abundance, the 
dashed lines ±0.2 dex of the adopted abundance. The CS 22891-171 spec- 
trum has a S/N per resolution element of 58, about average for our spectra. 

based on the differences in model atmospheres. 

In order to test this, and to show that our [Fe/H] values 
derived from Ca II K are good to within about 0.15 dex, we 
studied a subsample of the stars that had equivalent widths 
(EWs) based on high-resolution spectra available in the liter- 
ature. These EWs were run through our analysis, adopting 
our model atmosphere parameters; the results are shown in 
Table 5. The [Fe/H] derived from the Call K line, with an as- 
sumed [Ca/Fe] of +0.3, is in good agreement with the [Fe/H] 
derived from the EWs of Fe I lines, except for CS 22 1 83-03 1 , 
whose [Fe/H] in our analysis is 0.4 dex higher than in the EW 
analysis. The abundance ratios of [Ca/Fe], [Mg/Fe], [Ti/Fe] 
and [Fe II/Fe I] are also quite reasonable, although our higher 
adopted gravity for CS 22947-187 means that we derive much 
higher abundances for the ionized species than for the neutral 
species. From Table 5, we conclude that basing our model at- 
mosphere metallicities on the Ca II K line results in an offset 
of 0.12 dex with a rms scatter of ~ 0.16 dex, which is in line 
with the expected uncertainty due to our fits of the Ca II K re- 
gion. Our final model atmosphere parameters are summarized 
in Table 3. 

3.3. C and N Measurements and Upper Limits 

We measured C and N abundances by creating a synthetic 
spectrum for each star and comparing it with the observed 
spectra (Figure 6). When the NH and CH regions were syn- 
thesized, the abundance ratios of [Mg/Fe], [Ca/Fe] and [Ti/Fe] 
were set to be +0.3, +0.3, and +0.2, based on the Cayrel 
et al. (2004) results for the enhancements in the a-elements 
for metal-poor halo stars. Other elemental abundance ratios 
may also differ from solar in these stars, but their lines are ei- 
ther not strong or sufficiently plentiful to noticeably affect the 
synthesis. We calculated the molecular equilibrium among 
the species H2, CH, NH, OH, C2, CN, CO, N2, NO, O2, H2O, 
and CO2. Therefore, our derived C and N abundances could 
depend on the adopted O abundance. We set the [O/Fe] at 
+0.4 (Fulbright & Kraft 1999) in our syntheses, however the C 
and N abundances are not sensitive to the adopted [O/Fe] un- 
less it exceeds +1.5, and then only at the 0.01-0.02 dex level. 
The synthesis of the NH band does not depend on the C abun- 
dance, which was one of the advantages of this program. The 
random uncertainty from continuum placement, S/N, and im- 
perfect line lists is estimated to be 0.2 dex for the C, N and 
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TABLE 4 

Model Atmosphere Parameter Comparison 



Star 


Tgff 


log g 


[m/H] 


£, 


Source 


BD-18° 5550 


4750 


1.4 


-3.06 


1.80 


Cayrel et al. (2004) 




4806 


1.72 


-2.89 


1.91 


This study 


HD 122563 


4600 


1.10 


-2.82 


2.00 


Cayrel et al. (2004) 




4500 


1.30 


-2.74 


2.5 


Westin et al. (2000) 




4615 


1.27 


-2.47 


2.10 


This study 


HD 160617 


5967 


3.79 


-1.77 


1.50 


Jonsell et al. (2005) 




5931 


3.77 


-1.79 


1.50 


Akerman et al. (2004) 




5999 


3.74 


-1.36 


1.31 


Gratton et al. (2000) 




5882 


3.69 


-1.69 


1.45 


This study 


HD 186478 


4700 


1.30 


-2.59 


2.00 


Cayrel et al. (2004) 




4831 


1.78 


-2.63 


1.89 


This study 


CS 22183-031 


5270 


2.8 


-2.93 


1.20 


Honda et al. (2004b) 




5196 


2.76 


-2.79 


1.59 


This study 


CS 22892-052 


4850 


1.50 


-2.97 


2.50 


Norrisetal. (1997) 




4790 


1.60 


-2.92 


1.80 


Honda etal. (2004b) 




4850 


1.60 


-3.03 


1.90 


Cayrel et al. (2004) 




4760 


1.30 


-3.10 


2.29 


McWilliam et al. (1995a) 




4861 


1.86 


-2.76 


1.86 


This study 


CS 22947-187 


5160 


1.30 


-2.6 


2.26 


McWilliam et al. (1995a) 




5489 


3.44 


-2.25 


1.47 


This study 


CS 22950-046 


4640 


0.85 


-3.5 


2.68 


McWilUam etal. (1995a) 




4730 


1.30 


-3.30 


2.02 


Carretta et al. (2002) 




4604 


1.25 


-3.29 


2.11 


This study 


CS 22958-042 


6250 


3.50 


-2.85 


1.50 


Sivrani et al. (2006) 




6224 


3.86 


-2.65 


1.44 


This study 


CS 22968-014 


4840 


1.80 


-3.5 


1.90 


McWilliam et al. (1995a) 




4850 


1.70 


-3.56 


1.90 


Cayrel et al. (2004) 




4892 


1.93 


-3.30 


1.83 


This study 


CS 29497-030 


6650 


3.50 


-2.80 


2.00 


Sivarani et al. (2004) 




7000 


4.10 


-2.57 


1.90 


Ivans et al. (2005) 




7163 


4.20 


-2.20 


1.43 


This study 


CS 30314-067 


4400 


0.70 


-2.85 


2.50 


Aold et al. (2002d) 




4476 


0.96 


-2.67 


2.25 


This Study 


CS 31085-090 


6050 


4.50 


-2.85 


1.00 


Sivrani et al. (2006) 




5972 


3.73 


-2.75 


1.44 


This study 



TABLE 5 

Metallicity Comparison 



Our Synthesis Abundances Derived with Literature EWs and Our Model Atmospheres 

Star [Call/H] [Fe/H] [Cal/H] [Fel/H] [Cal/Fel] [Mgl/Fel] [Til/Fel] [Till/Fel] [Fell/Fel] EW Source 



CS 22183-031 


-2.49 


-2.79 


-2.76 


-3.19 


0.43 


0.72 


0.68 


0.58 


+0.13 


1 


CS 22947-187 


-1.95 


-2.25 


-1.77 


-2.11 


0.34 


0.44 


0.30 


0.96 


+0.59 


2 


CS 22892-052 


-2.46 


-2.76 


-2.65 


-2.94 


0.30 


0.46 


0.20 


0.16 


+0.00 


3 


CS 22950-046 


-2.99 


-3.29 


-3.23 


-3.46 


0.23 


0.63 


0.19 


0.22 


+0.00 


4 


CS 22968-014 


-3.00 


-3.30 


-3.35 


-3.38 


0.03 


0.47 


0.10 


0.06 


-0.05 


5 


CS 29497-030 


-1.90 


-2.20 


-1.85 


-2.30 


0.45 


0.66 


0.65 


0.47 


-0.15 


6 


CS 30314-067 


-2.37 


-2.67 


-2.61 


-2.71 


0.10 


0.69 


0.22 


0.42 


-0.06 
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(2005), (7)Aoki et al. (2002c) 



Ca syntheses based on by-eye fits to the data with different 
abundances and continuum placements. 

Our resolution is sufficiently low that a determination of a 
lower limit on N for some of our program stars relies on the 
lowest point in the spectra in the region of the NH feature. We 
found a lower limit by selecting the lowest pixel value within 
a window 0.8 A wide, centered on the expected position of 
the maximum N absorption. The rms scatter was determined 
from the S/N, and a 4-cr offset downward was added to that 
data point. Then, the lowest points in many NH syntheses 
were found, and interpolated in order to yield a lower limit on 
the N abundance. 



3 .4. Three -Dimensional and Non-LTE Effects 

Three-dimensonial effects are likely to be very important in 
deriving correct C and N abundances. Preliminary calcula- 
tions in red giant atmospheres for the NH lines at 3360 A and 
the CH lines at 4305 A indicate that the corrections can be 
as large as 0.4 and 0.5 dex, respectively (Garcia-Perez, priv. 
comm.) and 0.8 dex and larger for the most iron-poor stars 
(e.g. Collet et al. 2006). Because we are primarily concerned 
with [C/N] ratios in this paper, our results are not as sensi- 
tive to these corrections, and we ignore them for now. If 3-D 
corrections of this magnitude are appUed to our [C/Fe] val- 
ues, in many cases the values will drop below [C/Fe] = 0, 
and one might wonder if we still have a useful sample. How- 
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ever, in that case, [C/Fe] for metal-poor field stars in general 
will also drop by the same amount, and our stars will still be 
enriched relative to normal stars. Hence the assumption that 
C-enrichment is due to AGB stars still would stand. 

3.5. Error Analysis 

There are two major sources of error that we need to quan- 
tify. The first is the choice of molecular parameters, including 
the absolute oscillator strengths and dissociation potentials, 
which were discussed in § 3.1. The second is the choice of 
model atmosphere parameters, in particular Teff. Our sam- 
ples includes stars with a wide range of effective temperature 
and gravity, therefore we selected three stars covering these 
ranges, CS 22887-042, CS 22947-187, and CS 29493-090, 
for the error analysis. The changes in loge for C, N and Ca 
are listed in Table 6. The dependence on and metallicity 
are much smaller than those for Teff and log g, and hence will 
be ignored. The increase in log e with increasing temperature 
is offset by a decrease with increasing gravity when the ac- 
tual model atmosphere parameters are used. To account for 
the expected correlations between error sources, we adopt the 
following equation: 



'toga ■ 



where asyn is the abundance error associated with the synthe- 
sis and is 0.2 dex. In this case, where log g is derived from 
Teff, we write trriogg as 



(2) 



/ (91ogg 
I dT 



(3) 



The slope of the log g-Tea relation depends on the evolu- 
tionary state of the star We found slopes of 2.33 x 10"^ for the 
giant CS 29493-090, 1.27 x 10"^ for the subgiant CS 22947- 
187, and 9.33 x lO"'* for the turnoff star CS 29493-090. The 
errors in [X/H] were calculated using equation 2, while the 
the abundance ratio errors for [C/N], [C/Fe] and [N/Fe] were 
calculated using a modification of Equations A19 and A20 
from McWilliam et al. (1995a). Our uncertainties for turnoff, 
subgiant, and giant stars are summarized in Table 7. 

4. RESULTS 

The abundances of N and C, derived from the NH, CN, and 
CH lines, are listed in Table 8 and plotted in Figure 7. Twenty- 
one of the stars in our sample are C-mild to C-rich and are per- 
tinent to the problem of the "missing" NEMP stars. The seven 
often-observed calibration stars (CS 22892-052, CS 22968- 
014, CS 22950-046 and the HD/BD stars) are not shown in 
Figure 7. It is clear from the low Teff of CS 22950-046 that 
its low [C/N] ratio is due to internal mixing (Gratton et al. 
2000; Spite et al. 2005), and it is similar to the other low Teff 
calibration giants, HD 122563 and HD 186478. 

Two facts are immediately apparent from inspection of Fig- 
ure 7.(1) We have mild C-enhancements for many stars in our 
sample, and have managed to fill in the +0.5 < [C/Fe] < +1 .0 
region of the diagram and (2) our [C/N] values are not con- 
sistent with those expected to be produced by HBB in AGB 
stars, and indeed are rather similar to previous results. There- 
fore, it appears that an observational bias is not the solution to 
the scarcity of NEMP stars. 




HE 1337-3326 



A-HE 1031-0020 



CS 29528-041 



■ This study <r 
ft s-process 
A no s— process 



CS 23949-037 



CS 30322-033 



[C/Fe] 

Fig. 7. — [C/N] vs. [C/Fe] for our sample of stars (filled squares) as well as 
literature stars with (stars) and without (triangles) s-process enhancements. 
The lines are the same as in Figure 2. The [C/N] ratios for the most iron-poor 
stars are indicated by rightward pointing arrows because their [C/Fe] ratios 
(HE 0107-5240, [C/Fe] = +3.71, and HE 1327-2326, [C/Fe] = +4.26) are 
off the scale of the plot. 

Table 9 lists our derived [C/Fe] and [C/N] ratios compared 
with available literature values for stars in common, while 
Table 10 lists the additional literature values used in Fig- 
ure 7. There are three stars in the non-calibration sample 
with [C /Fe] < 0. However, examination of Table 9 shows 
that our [C/Fe] ratios are in general lower than the literature 
values. Because these stars were chosen out of a large sam- 
ple with consistent C-abundance determinations (Rossi et al. 
2005), we believe they are still C-enhanced relative to the field 
population. 

There are four stars with literature [C/N] values in the in- 
teresting low range of [C/N] in Figure 7; these deserve some 
special discussion. Based on the observations in Figure 7, 
we define the class of NEIVIP stars as having [N/Fe]>0.5 
and [C/N]< -0.5. This marks stars that have distinctly lower 
[C/N] ratios than the rest of the CEIVIP stars, although we note 
that to agree with the predictions of HBB, these stars would 
need to have [C/N] ~ -1. Two of these, CS 29528-041 and 
HE 1031-0020 have properties similar to the NE]V1P stars we 
were expecting to find in much greater numbers. Of the other 
two, CS 22949-037 is not s-process rich (but it is oxygen rich), 
and CS 30322-023 is probably an intrinsic thermally pulsing 
AGB star (Masseron et al. 2006). The statistics for the pres- 
ence of NEIVIP stars in the literature confirm the rarity that led 
to the lack of them in our sample. 

• CS 29528-041: Sivrani et al. (2006) report on the ele- 
mental abundances of this C-rich, but even more N-rich 
star. It has the expected abundance ratios for NEIVIP 
stars: [Fe/H] = -3.3, [C/Fe] = +1.6, [N/Fe] = +3.0, 
hence [C/N] = -1 .4. Interestingly, this star also has de- 
tectable lithium, albeit at a value well below the Spite 
Plateau, A(Li) =1.7. It is s-process-element rich as 
well, with [Ba/Fe] = +1.0, which indicates that more 
massive AGB stars can produce the s-process. 

• HE 1031-0020: Cohen et al. (2006) found a [C/N] 
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TABLE 6 

A Abundance for C, N, and Ca with Model Atmosphere Changes 



Star 


Aloge(C) 


Aloj 


;e(N) 


Aloge(Ca) 




A Teff Alog g 


A Teff 


Alogg 


A Teff Alog g 




+150 K +0.3 dex 


+150 K 


+0.3 dex 


+150 K +0.3 dex 


CS 22887-048 


0.250 -0.125 


0.325 


-0.125 


0.125 0.000 


CS 22947-187 


0.275 -0.100 


0.275 


-0.125 


0.015 -0.050 


CS 29493-090 


0.310 -0.125 


0.300 


-0.225 


0.100 -0.025 



TABLE 7 
Summary of Errors 



Star 


Teff range 


^rClogeCO) 


<T(loge(N)) 


(T(loge(Fe)) 


<7([C/Fe]) 


(7([N/Fe]) 


o-([C/N]) 


turnoff 


Teff> 6400K 


0.30 


0.36 


0.24 


0.31 


0.35 


0.31 


subgiant 


5400K<Tetf< 6400K 


0.30 


0.31 


0.24 


0.30 


0.32 


0.30 


giant 


Teft< 5400K 


0.25 


0.25 


0.21 


0.29 


0.25 


0.41 



of —0.85 for this s-process-element rich star, which is 
one of the less C-enhanced stars in their sample, with 
[C/Fe]=+1.6. 

CS 22949-037: Norris et al. (2002) showed that this 
star has a [C/N] ratio of -1.25, while Depagne et al. 
(2002) found an extremely high [O/Fe] ratio of +1.97. 
Models of EMP or Z = low-mass AGB stars based 
on standard mixing assumptions (no convective over- 
shooting at bottom of He-shell flash convection zone) 
do predict a range of O overabundances depending on 
mass (Siess et al. 2002; Herwig 2004b). ForZ = 0.0001, 
[O/Fe] can be up to 1.7 dex for low-mass models. Like 
C and N, the O production in AGB stars is primary, and 
therefore larger overabundances are expected with low 
metal content. However, the large N overabundance ob- 
served in this star is only predicted for more massive 
AGB stars for which the O overabundance is proba- 
bly smaller. This, and the lack of enrichment in the 
s-process elements, supports the hypernova hypothesis 
for this star. 

CS 30322-023: Masseron et al. (2006) have argued that 
this extremely low metallicity star, with [Fe/H] = -3.5, 
and with a clear s-process-element signature, is likely 
to be an example of an intrinsic AGB star, caught dur- 
ing the brief stage of evolution when thermal pulses 
are occurring. This star is mildly carbon-rich, with 
[C/Fe] = +0.6, while the nitrogen abundance ratio is 
quite high, [N/Fe] = +2.8. We would not classify this 
an NEMP star, however, if it is an intrinsic AGB star, 
because there is no binary mass transfer, in analogy to 
the CEMP-s stars. Note, however, that if the identifica- 
tion of this star as a TP- AGB is correct, its mass cannot 
be high, but rather, it must be on the order of 0.8 Mq. 
Since at masses this low, HBB is not expected to oc- 
cur, Masseron et al. (2006) suggest an unknown mixing 
process that happens in VMP stars is responsible. 



DISCUSSION: WHY THE LACK OF OBSERVED 
NEMP STARS? 



We expected to find that between 12% and 35% of our 
sample would turn out to be NEMP stars (see § 5.2), depend- 
ing on our assumption concerning the mass ratios in the binary 
progenitors. Instead, we found no candidates in our sample of 
21 stars, and only two in the recent literature. The presence 
of these two [C/NJ ^ -1.00, s-process-rich stars suggest that 
NEMP stars do in fact exist, and their abundance patterns are 
as expected from the models, just not in the predicted numbers 
that should be found. However, there remain several outstand- 
ing questions concerning the nature of our sample which we 
should examine before asking if models of primary nitrogen 
production by intermediate-mass AGB stars are incorrect. 

5.1. Are we observing stars that have been polluted by AGB 

stars? 

AGB stars are by far the most likely culprits for the pollu- 
tion of our sample. The majority of CEMP stars are CEMP-s 
stars, and much observational evidence, reviewed in the In- 
troduction, points to their being the result of AGB mass trans- 
fer. However, it is possible that stars with lower [C/Fe] values 
have a smaller percentage that are s-process-rich than those 
with [C/Fe] > 1.0. We have information on the s-process 
abundances for four stars in the sample: CS 22947-187 and 
CS 31085-090 are s-process-rich, while CS 30314-067 and 
CS 22958-042 are not and there appears to be no correlation 
with C-richness. We are currently obtaining follow-up high- 
resolution spectra for many of stars in our sample to check for 
radial-velocity variations and s-process enhancements. 

5.2. Are we observing stars that have been polluted by 
intermediate-mass AGB stars? 

Although we are confident that most of our stars have been 
polluted by AGB stars, we have considerably less information 
on the distribution of masses of those AGB stars. Our initial 
estimate was based on the mass ratio (q) distribution deduced 
by Pinsonneault & Stanek (2006) from a variety of studies, 
particularly the data on massive binaries in the Small Mag- 
ellanic Cloud from Harries et al. (2003) and Hilditch et al. 
(2005). Their most Ukely distribution of q fractions is that 
45 % of binary systems are "twins", with the primary and 
secondary masses the same to within 5%, while the q of 
the other 55 % of systems are matched by a flat mass dis- 
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TABLE 8 
C, N, AND [Fe/H] for Our Sample 



Star loge(C) [C/Fe] loge(N) [N/Fe] [Fe/H] [C/N] 
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6 
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40, 


.25 
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.62 
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.96 


-2, 


.39 


> 


-0 


.71 


CS 22949-008 


7, 


.79 


+1 


.15 


6, 


.43 




+0, 


.30 


-1, 


.92 




+0 


.85 


CS 22950-046 


4, 


.77 


-0 


.50 


6, 


.16 




+0, 


.30 


-3, 


.29 




+0 


.80 


CS 22958-042 


8, 


.31 


+2, 


.40 


<7, 


.29 
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-t-1, 


.89 


-2, 


,65 


> 


-fO 


.51 


CS 22958-083 


6 


.41 


+0 


.35 


<6, 


.51 


< 


40, 


.96 


-2, 


,50 


> 


-0 


.61 


CS 22960-053 


6 


.63 


+1 


.15 


6. 


.12 




+1, 


.15 


-3. 


.08 




+0 


.00 


CS 22968-014 


5, 


.36 


+0, 


.10 


<4, 


.88 


< 


+0, 


.13 


-3, 


.30 


> 


-0, 


.03 


CS 29493-090 


5, 


.99 


+0, 


.25 


5, 


.68 




+0, 


.45 


-2, 


.82 




-0, 


.20 


CS 29495-042 


7, 


.26 


+1, 


.00 


6, 


.25 




+0, 


.50 


-2, 


.30 




+0, 


.50 


CS 29497-030 


8, 


.36 


+2, 


.00 


<7, 


.99 
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+2, 


.14 


-2, 


.20 


> 


-0, 


.14 


CS 29512-073 


7, 


.51 


+1, 


.05 


5, 


.95 




+0, 


.00 


-2, 


.10 




+1, 


.05 


CS 30314-067 


6, 


.14 


+0, 


.25 


5, 


.88 




+0, 


.50 


-2, 


.67 




-0, 


.25 


CS 31062-041 


6 


.51 


+0, 


.25 


6 


,15 




+0, 


,40 


-2, 


,30 




+0, 


.15 


CS 31080-095 


7, 


.91 


+2, 


.10 
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.26 
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TABLE 9 

Comparison with Literature Values 



Star 


[C/H] 


[N/H] 


[C/N] 


[C/Fe] 


Source 


BD-18°5550 


-3.08 


-3.02 




-0.02 


Spite et al. (2005) 




-3.09 






-0.20 


This study 


HD 122563 


-3.29 


-1.72 


-1.57 


-0.47 


Spite et al. (2005) 




-3.20 


-1.65 


-1.55 


-0.46 


Westin et al. (2000) 




-3.42 


-1.47 


-1.95 


-0.95 


This study 


HD 160617 


-1.82 






+0.03 


Akerman et al. (2004) 




< -1.69 


-0.34 


<-1.35 


<+0.30 


Lau-d (1985) 




-2.29 


-1.29 


-1.00 


-0.60 


This study 


HD 186478 


-2.89 


-1.57 


-1.32 


-0.30 


Spite et al. (2005) 




-2.83 


-1.93 


-1.00 


-0.20 


This study 


CS 22892-052 


-2.06 






+0.98 


McWilliam et al. (1995a) 




-2.26 


-2.22 


-0.04 


+0.88 


Sneden et al. (2003) 




-2.16 


-2.91 


+0.75 


+0.60 


This study 


CS 22958-042 


+0.14 


-1.00 


+1.14 


+3.01 


Sivrani et al. (2006) 




-0.25 


< -0.76 


>+0.51 


+2.40 


This study 


CS 30314-067 


-2.4 


-1.7 


-0.70 


+0.50 


Aoki et al. (2002c) 




-2.42 


-2.17 


-0.25 


+0.25 


This study 


CS 31085-095 


-0.31 


-2.40 


+2.09 


+2.56 


Sivrani et al. (2006) 




-0.65 


-1.49 


>+0.84 


+2.10 


This study 



tribution. We assumed that the observed stars, being main- 
sequence tumoff, subgiant, giant or horizontal branch stars, 
all have main-sequence masses close to 0.9 M©. Therefore, 
their "twin" AGB stars would be < 1 Mq and would have 
created high [C/N] ratios. Therefore, the NEMP stars must 
come from the flat mass distribution. We adopt 3.5 Mq as the 
dividing line between those AGB stars that undergo HBB and 
those that do not. FinaUy, the maximum mass of an AGB star 
at [Fe/H] = -2.3 is 7.5 Mq (Poelarends et al., in prep). Com- 
bining all of this information, we find that there should be 
roughly 33 % NEIVIP stars and 67 % CElVIP-s stars contained 



in our sample. 

There may be difficulties with our predictions. For exam- 
ple, while Pinsonneault & Stanek (2006) consider the mass 
ratios using studies that cover a wide range of masses, none 
of the studies is able to cover a wide range of q. The SMC bi- 
nary studies only reach a mass ratio of ~ 0.4, which is larger 
than minimum q we are interested in, 0.12, and it is conceiv- 
able that the distribution will change as the mass difference 
between the two components becomes more extreme. Such 
studies are difficult because of the different luminosities and 
timescales of the two kinds of stars, but we do have some 
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TABLE 10 
Additional Literature Values 



Star loge(C) loge(N) [Fe/H] [C/N] Source s-process? 



Original Literature Sample 



moT~7 ^'\f\^ 


6.86 


5.35 


—2.72 


+1.00 


\ ^1,; „1 /OAAO^\ 

AoKi et al. (2002c) 


no 


CS 22560-074 


7.93 


6.02 


—1.93 


+1.40 


Aoki et al. (2002d) 


yes 


Co 22ool-U3o 


8.46 


6.99 


—2.06 


+0.96 


Preston & Sneden (2001) 


yes 


Co zioyo-Ui/ 


8.51 


6.70 


—2.25 


+1.30 


AoKi et al. (2002d) 


yes 


Co z2y4z-uiy 


7.92 


6.01 


-2.64 


+1.70 


Aoki et al. (2002d) 


yes 


/^c oor»/io /^^OT 

Co 22940-02/ 


7.96 


7.25 


-2.57 


+0.20 


Aoki et al. (2()02c) 


no 


CS 22957-027 


7.80 


6.45 


-3.38 


+0.84 


Aoki et al. (2()02b) 


no 


/^c or\/im /T) /I 

Co 29497-034 


7.56 


7.45 


-2.90 


-0.40 


Hill et al. (2000) 


yes 


Co 29490-043 


6.70 


6.50 


-3.75 


-0.31 


Aoki et al. (2()02b) 


no 


/^c onc/^^o r\r\o 

Co 29502-092 


6.76 


5.95 


-2.76 


+0.30 


Aoki et al. (2()02c) 


no 


/^c or\coz; i i a 

CS 29526-1 10 


8.38 


7.07 


-2.38 


+0.80 


Aoki et al. (2002d) 


yes 


CS 30301-015 


7.52 


6.01 


-2.64 


+1.00 


Aoki et al. (2()02d) 


yes 


CS 31062-012 


8.11 


6.70 


-2.55 


+0.90 


Aoki et al. (2002d) 


yes 


CS 31062-050 


8.24 


6.93 


-2.32 


+0.80 


Aoki et al. (2002d) 


yes 


HD 196944 


7.51 


7.10 


-2.25 


-0.10 


Aoki et al. (2002d) 


yes 


JrLh UU24-2523 


8.44 


7.43 


-2.72 


+0.50 


LucateUo et al. (2003) 


yes 


LP 625-44 


8.00 


6.30 


-2.71 


+1.19 


Aoki et al. (2001) 


yes 




7.96 


7.03 


-2.74 


+0.42 


Aokietal. (2001) 


yes 






New and Updated 


Literature Sample 


T TT~^ /■^^'^ 1 ^ 1/1/11 

Hb 0012-1441 


7.66 


6.05 


-2.52 


+1.10 


Cohen et al. (2006) 


yes 


I TT7 A A C O AO /I A 

Hh 0050-0244 


7.76 


6.95 


-2.75 


+0.30 


Cohen et al. (2006) 


yes 


Hh 0107-5240 


6.81 


5.22 


-5.46 


+1.08 


Christlieb et al. (2004) 


no 


I TT? /^M /ll /^^ /I /I 1 

Hh 0143-0441 


8.26 


7.35 


-2.31 


+0.40 


Cohen et al. (2006) 


yes 


HF 021 2-05^7 


8.06 


6.75 


-2.27 


+0.80 


Cohen et al. (2006) 




HE 0336+0113 


8.16 


6.85 


-2.68 


+0.80 


Cohen et al. (2006) 


yes 


HE 1 150-0428 


7.66 


7.15 


-3.30 


+0.00 


Cohen et al. (2006) 


no 


HE 1410+0213 


8.16 


6.85 


-2.16 


+0.10 


Cohen et al. (2006) 


no 


HE 1434-1442 


8.16 


6.95 


-2.39 


+0.70 


Cohen et al. (2006) 


yes 


HE 1509-0806 


7.66 


7.25 


-2.91 


-0.10 


Cohen et al. (2006) 


yes 


HE 2158-0348 


7.76 


6.75 


-2.70 


+0.50 


Cohen et al. (2006) 


yes 


HE 2232-0603 


7.96 


6.55 


-1.85 


+0.90 


Cohen et al. (2006) 


yes 


HE 2356-0410 


7.66 


6.75 


-3.07 


+0.40 


Cohen et al. (2006) 


no 


HE 1327-2326 


6.99 


6.83 


-5.73 


-0.37 


Aoki et al. (2006) 


no 








Nitrogen-Rich Stars 




CS 29529-041 


6.70 


7.57 


-3.32 


-1.38 


Sivrani et al. (2006) 


yes 


CS 22949-037 


5.82 


6.96 


-3.79 


-1.65 


Norris et al. (2002) 


no 


CS 30322-023 


5.60 


7.20 


-3.39 


-2.11 


Masseron et al. (2006) 


yes 


HE 1031-0020 


7.36 


7.55 


-2.86 


-0.70 


Cohen et al. (2006) 


yes 



information on the q distribution of binaries with solar-type 
primaries (Duquennoy & Mayor 1991). For the longer period 
(P > 3000 days) binaries, the data suggested a secondary star 
drawn from an IMF-weighted distribution. However, for the 
shorter period binaries, Mazeh et al. (1992) found a distribu- 
tion that was much flatter, with perhaps a rise toward q= \. 
For lack of other information, if we make the assumption that 
the mass of the companion of the observed stars is drawn from 
a Salpeter (1955) mass function between IMq and 7.5 M0, 
we find that there should be 1 NEMP for every 6.9 CEMP-s 
stars. This is a much smaller ratio, but still not small enough 
to be supported by the observational data. All of these num- 
bers were obtained assuming that there was not a bias towards 
finding NEMP stars in the C-mild sample; if such a bias ex- 
isted, the disagreement with theory would only be strength- 
ened. 

5.3. Could binarity be affecting our results? 

The fact that these AGB stars were in binary systems may 
have affected their nucleosythesis yields by increasing their 
mass-loss rates and ending the production of N by HBB. This 
could be due to tidal synchronization between rotational and 
orbital motions, as observationally demonstrated for RGB 
stars by De Medeiros et al. (2002). Even after HBB is ter- 
minated, dredge-up can continue if the minimum envelope 



mass for dredge-up is less than that for HBB, and the result- 
ing [C/N] ratio may be altered. Frost et al. (1998) argued that 
the presence of high- luminosity C stars in the Large Magel- 
lanic Cloud could be explained by continued dredge-up after 
HBB has ceased. It is not clear if this same reasoning would 
apply to our present sample, which has much lower metallic- 
ity. However, if tidal synchronization leads to a spin-up mix- 
ing, N production should be enhanced, possibly balancing the 
effect of enhanced mass loss. Without detailed models we 
can only speculate, but it is not immediately clear that bina- 
rity could inhibit N production at the bottom of more massive 
AGB stars. 

5.4. Could mass ratios affect mass transfer? 

Jorissen & Boffin (1992) summarize the situation regarding 
mass transfer in both Roche-lobe overflow and wind models. 
For large mass ratios, only wind models were stable, because 
the large mass-transfer rates in overflow scenarios puff up the 
secondary star and lead to a contact binary situation (Tout & 
Hall 1991). Vanture ( 1 992) argued that the mass transfer from 
a wind would not be sufficient to create CH stars, and there- 
fore only in binary systems with Roche-lobe overflow would 
the nucleosynthesis of the AGB star be preserved in the still- 
shining secondary. Much work remains to be done in under- 
standing mass-transfer scenarios, but if NEMP stars can only 
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be formed by wind accretion, then their production could be 
far less efficient than for CEMP stars. 

5.5. Could mass transfer rates be affected C- vs. 

N- richness? 

Our understanding of mass loss in AGB stars is seriously 
incomplete. However, one variable that may be important is 
the dust-to-gas ratio, with mass-loss rates increasing for in- 
creasing amounts of dust. Carbon is a major contributor to 
dust in AGB stars, and the dust-to-gas ratio has been param- 
eterized by the C-richness (e.g., Arndt et al. 1997). There- 
fore, since HBB stars produce less C-rich AGB envelopes, 
their mass loss may be lower and therefore the number of bi- 
nary companions that are enriched enough to be included in 
our sample may be small. However, a lower mass-loss rate, 
given enough time, could still conceivably add up to the same 
amount of enrichment. Additionally, Arndt et al. (1997) found 
only a small dependency on C-enhancement. 

5.6. Does hot bottom-burning occur in intermediate mass 
AGB stars? 

Another explanation for the absence of NEMP stars in our 
sample could be the inability of more massive AGB stars 
(maybe only at very low metallicity) to produce N via HBB. 
The efficiency of envelope convection, parameterized by the 
mixing length parameter, determines the efficiency of HBB 
in stellar evolution models. A larger convective efficiency 
would imply even more efficient HBB, but could decrease 
the dredge-up efficiency, in particular if significant convective 
overshoot is present. This would limit the ability of massive 
AGB stars to produce primary N. A smaller convective effi- 
ciency would lead to less efficient HBB, and a smaller mass 
range for efficient N production, lowering the expected num- 
ber of NEMP stars. However, we have not performed calcu- 
lations with a range of convective efficiencies, and therefore 
we cannot quantify these possible effects. Stellar evolution 
modeling of HBB in massive AGB stars also predicts the pro- 
duction of Li through the Cameron & Fowler (1971) mecha- 
nism. This property of HBB is observationally confirmed by 
the discovery that almost all of the luminous AGB stars with 
C/O < 1 in the Magellanic Clouds are lihtium rich (e.g. Smith 
et al. 1995). 

5.7. Could extra-mixing be affecting our results? 

The main goal of our study was to establish whether there 
existed a possible observational selection bias against the dis- 
covery of NEMP stars. However, we also noted that previ- 
ous literature data, as well as our new objects have C/N ratios 
that are too low compared to standard low-mass AGB models 



without HBB. Our 2 M© and 3 Mq model sequences produce 
large amounts of C, but no N. The most likely solution to this 
discrepancy would be an additional mixing process, such as 
cool-bottom processing, which has been studied in great detail 
in the context of abundance anomalies in globular cluster stars 
(Wasserburg et al. 1995; Denissenkov & VandenBerg 2003). 
Such an extra-mixing process below the bottom of the con- 
vective envelope could also operate in low-mass AGB stars, 
leading to mild CN cycling that produces some N, but not as 
complete as in HBB (NoUett et al. 2003). We beUeve that the 
present observational data support such a scenario. We plan to 
test this idea quantitatively in the future with stellar evolution 
models that include this process. 

5.8. Final Thoughts 

It is probable that the formation of C-rich stars involves a 
delicate balancing act in which the N-rich stars are disadvan- 
taged. The most likely possibiHties are the binary mass ra- 
tios in the early Universe or the process of mass transfer for 
large mass-ratio systems. With ever larger samples of metal- 
poor stars with sufficient spectroscopic data becoming avail- 
able, we expect to continue studies of the relative fractions of 
CEMP vs. NEMP stars. To confirm that the majority of the 
CEMP stars in our sample were indeed polluted by AGB pro- 
genitors, we and other authors are obtaining high-resolution 
spectra for many of these stars to look for s-process enhance- 
ments, and to check for detectable radial-velocity variations. 
In addition, other elements, such as Li, Na, F, and the isotopes 
of Mg, which are affected by nucleosynthesis in AGB stars, 
will be examined in due course. 
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